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At the nucleus of a superconducting vortex is a small, circular core region where superconductivity 
is destroyed. Like an atomic nucleus, this core may become deformed, and such distortions can 
have important consequences in non-equilibrium situations. Using Ginzburg-Landau theory, we 
have investigated this phenomenon for vortices in the presence of artificial defects. We show that 
when a vortex approaches the vicinity of a defect, an abrupt transition occurs in which the vortex 
core develops a "string" extending to the defect boundary, while simultaneously the supercurrents 
and associated magnetic flux spread out and engulf the defect. The energetics of stretching the 
string determines the pinning behavior of the vortex. Experimental consequences of these strings 
are discussed. 



PACS numbers; 74.20.-z, 74.60. Ge, 74.76.-w, 68.55.Ln, 61.46.-fw 



One of the most important properties of superconduc- 
tors is their ability to carry currents without dissipation, 
allowing them to generate large magnetic fields. Many 
superconductors allow fields to penetrate in bundles of 
quantized magnetic flux, with associated whirlpools of 
current known as vortices. When these vortices are mo- 
bile, they spoil the perfect conductivity that make su- 
perconductors so useful. The quest to increase the max- 
imum dissipationless current Jc that a superconductor 
may carry has thus fueled intense study of vortex pin- 
ning. 

In recent years, pinning environments of artificially 
fabricated nanoscale defect arrays have been developed in 
hopes of better understanding and improving the pinning 
properties of superconductors. Some of the earlier con- 
tributions have involved macroscopic measurements (e.g. 
Jc, magnetization) performed on periodic "antidot" ar- 
rays [^]-|3| . The antidot regions contain material which is 
rendered non-superconducting. Pinning behavior may be 
studied in these periodic systems using various imaging 
techniques 

Much of the theoretical work on vortex pinning has 
employed numerical studies to focus |^ on the behavior 
of large collections of vortices under the influence of a 
driving force (supercurrent). Such studies usually em- 
ploy simplified pinning potentials, in part to make pos- 
sible simulations of large numbers of vortices, but also 
because information about pinning potentials at the mi- 
croscopic level is simply unavailable. A few studies PJlO| 
have focused on energy scales for pinning varying num- 
bers of fiux quanta to the defects, as well as defect- vortex 
potentials as derived from the London equation Jll| , p^ . 
However, the latter approach does not allow for varia- 
tion of the Cooper pair density, and in particular cannot 
correctly treat the vortex core. In the London approach, 
vortex cores are usually assumed to be rigid in shape, 
and interactions of vortices with their environments are 
determined by the core position as well as the distribu- 



tion of currents |ll| , |l2{ . Our results demonstrate that the 
vortex core in fact deforms dramatically near an artifi- 
cial defect: when the vortex center is sufficiently close to 
the defect, a string of suppressed order parameter devel- 
ops from the vortex position to the pinning center edge, 
while the currents and magnetic flux spread out over a 
large area (see Fig. 2 below). The pinning potential of 
the vortex turns out to be dominated by the string when 
it is present, as we discuss below. 

Methods- Our calculations focus on two dimensional 
arrays of artificial pinning centers in the form of holes 
in a bulk superconductor. Our goal is to find the low- 
est energy state of the system for specified locations of a 
superconducting vortex; from this we can construct a pin- 
ning potential. The appropriate description of the super- 
conducting state is in terms of Ginzburg-Landau theory, 
which focuses on a complex superconducting order pa- 
rameter ipi^^j for which lipl^ is proportional to the local 
density of superconducting electrons. Unlike the Lon- 
don theory, Ginzburg-Landau theory is valid at scales 
as small as the coherence length, ^. Written in terms 
of dimensionless variables, the Ginzburg-Landau energy 
functional is 



Egl = 



(1) 



In Eq. 1^, A is the vector potential, and the magnetic field 
B{r) = V X A. K = -I is the Ginzburg parameter, the 
ratio of the magnetic penetration depth A and the coher- 
ence length. In this work, we report on results obtained 
for K = 8. This choice of k is deep enough into the high 
K limit that, apart from scale factors in Eql and ip, the 
results vary little as k, is increased. 

To analyze the behavior of the vortex near the pinning 
center, we employ a mean field approach in which one 
minimizes Eql for a fixed vortex location to find V': A, 
and the current J. Our strategy for calculating -ip and 
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A self-consistently involves first holding ip fixed at some 
initial guess, and minimizing Egl with respect to A and 
B. Next, we fix A and B and minimize with respect to 
■0. These steps are iterated until changes in the variables 
become negligible. We implement this self-consistent ap- 
proach numerically by dividing the unit cell into a fine 
lattice of small unit cells. In this discrete scheme, V'(r) is 
replaced by ipij with ij specifying a grid point on a square 
lattice, while Aij and Jtj are defined on nearest-neighbor 
links between the grid points. Derivatives in Eq. |l] arc 
replaced by the corresponding finite differences. The re- 
sulting theory is very similar to lattice gauge theories 
studied in particle physics. To model the defect array 
as accurately as possible, one desires a fine grid; we find 
that with a 128 x 128 mesh our results are well-converged 
with respect to the discretization. 

To see how one minimizes Egl under the constraint 
of a specified vortex location, it is useful to write the 
current J in terms of the order parameter ip and the 
vector potential A. By minimizing Eql with respect to 
the vector potential and employing a Maxwell equation 
one has 
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In Eq. ^, we have used for the order parameter = 
|'0|e*'^. The familiar fluxoid quantization condition ||l^] 
arises from the requirement that the order parameter be 
single valued, i.e. ^ V0 • ds = 2Trny. Hence, in terms of 
J and A, 



for antidot systems the results obtained are quite similar 
to the large k results we report here. This may be un- 
derstood in terms of the effective magnetic penetration 
depth for a thin film, Ae/ / = — , which is typically much 
larger than the bulk value A [Q. This means that the 
energy stored in the magnetic field generated by the su- 
percurrents is quite small, so that the fact that the field 
varies as one moves out of the plane has little impact on 
the state of the system. The resulting energy functional is 
thus nearly identical to the bulk three dimensional case, 
with columnar antidots and vortices. 

Results- At large vortex-defect separations, the core 
has the usual compact structure with supercurrents lo- 
calized about it. Fig. ^ presents a perspective plot of |0|, 
as well as a vector plot of the currents. The distances 
shown are in units of the coherence length, ^. Our choice 
of a unit cell with side spanning 20 coherence lengths is 
typical of many of the nanoscale arrays studied exper- 
imentally. As the flux quantum nears the defect edge, 
it eventually reaches a critical distance dc (typically sev- 
eral ^, with precise value depending on the dot size and 
shape), where there is a sudden dramatic change. Fig. ^ 
illustrates the situation after the transition: the vortex 
core has developed a string extending from the flux quan- 
tum position to the defect edge; simultaneously the cur- 
rent now encircles the vortex-defect pair, and the mag- 
netic flux created by these currents spreads over a larger 
area. The string is energetically favorable because it al- 
lows the formerly dense current of the vortex to spread 
out (engulfing the defect in the process), thereby reduc- 
ing the kinetic energy of the state. 
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In the second part of Eq. ^, $5 is the total magnetic 
flux passing through the area of the contour, which we 
conveniently choose to be the small unit cell associated 
with the grid point ij, while n^iij) is the total number of 
"fluxoid quanta" contained in the contour of integration 
p^ . It is through n„ that the vortex location(s) in the 
full unit cell of the system may be fixed: riy = except 
at the grid points where we wish to place a vortex, for 
which Uy — 1. Armed with knowledge of |0| and some 
specified realization of n^{ij), one can solve for J and A 
via Eq. ^. Using the expression for the current given in 
Eq. ^, one obtains V0; inserting V0 and A into Eq. |l| 
yields an expression depending only on |0| and k, which 
we minimize with respect to |0| to obtain the order pa- 
rameter modulus. 

We note that the above method can be generalized to 
the case of a thin film superconductor. We have per- 
formed some calculations for such systems, and find that 
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FIG. 1. A perspective plot of the order parameter modulus, 
(a) and current image (b) just prior to the transition. 
Currents are localized about the vortex core, which has a 
compact structure. 



x/5 




5 10 15 

(b) 

FIG. 2. Order parameter perspective plot (a) and currents 
(b) just after the transition. The currents circulate about the 
vortex-defect pair, and the vortex core has a string extending 
all the way to the defect edge. 



When a net supercurrent flows across a superconduc- 
tor, a Lorentz force is exerted on vortices; if these move 
in response energy is dissipated. In an artificial defect 
array, the driving force can be balanced by a pinning 
force given by the gradient of Egl with respect to the 
vortex position. Hence, it is appropriate to regard Egl 
as the pinning potential. Fig. |^ shows this as a func- 
tion of distance from the unit cell edge, and it can be 
seen to have three distinct regions. For large separations, 
where the vortex core is compact, the pinning potential 
decreases relatively slowly. As the separation decreases 
it eventually crosses dc and the core-string structure ap- 
pears. Hysteresis in the calculations, indicated with ar- 
rows, suggests that the transition is first order. Below 
the transition one observes nearly linear behavior of the 
pinning potential, suggesting that the string carries an 
energy proportional to its length. The third region is 
announced by a discontinuous jump as the vortex is ab- 
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sorbed by the defect, followed by a perfectly flat region 
inside the defect. This abrupt jump is due to a sud- 
den transformation of the order parameter. Immediately 
prior to absorption, constant Cooper pair density con- 
tours near the flux quantum have a semicircular profile. 
As the vortex enters the defect, this "semi- vortex" van- 
ishes, removing a finite amount of energy for an infinites- 
imal change in the vortex position. 
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FIG. 3. Plot of the pinning potential showing broad linear 
regions and a gap in the potential at the defect boundary. The 
horizontal axis measures the distance of the vortex from the 
edge of the unit cell. Arrows and red circles indicate hysteresis 
in the calculation, a hallmark of a first order transition. 



Experimental implications- The profile of the pinning 
potential shown in Fig. ^ suggests a set of measurements 
one might perform to seek an experimental signature of 
the unique pinning phenomena discussed above. One 
possible test would involve an AC driving current of mag- 
nitude small enough to leave the vortices pinned by the 
strings. Such a current would allow the vortices to "rattle 
around" in the linear part of the pinning well illustrated 
in Fig. ^ and produce losses, whereas a DC current of 
equal magnitude would be dissipationless. Upon lower- 
ing of the temperature the vortices would be captured 
by defects, leading to lossless supercurrents for both AC 
and DC driving forces. An observation of these effects 
would yield indirect confirmation of the form of the pin- 
ning potential we find. Something like this may recently 
have been observed |l5). Another interesting possibil- 
ity is that a unique thermal depinning may occur as the 
temperature is increased in the regime of linear pinning. 
The presence of a string suggests that this may carry en- 
tropy at finite temperature, much as is the case of poly- 
mers. This entropy is proportional to the string length 
and temperature, and at high enough temperatures may 



overwhelm the energy per unit length found in our mean- 
field calculations. In analogy with polymer behavior [ p^ , 
this leads to unbounded growth of the string and effective 
depinning of the vortex. However, it is not clear whether 
the string remains sufficiently well-defined at the temper- 
atures necessary for proliferation that the polymer anal- 
ogy remains valid up to the transition. Further research 
into this possibility is currently underway. 

Ultimately the best confirmation of our results would 
involve imaging of the string. An interesting possibil- 
ity in this context is to fabricate arrays with two defects 
per unit cell in close proximity. When one magnetic flux 
quantum passes through each unit cell, we have found 
solutions in which a string develops between the defects. 
Fig. ^ is an image of j^/ij for an equilibrium configuration 
in a typical case. The existence of a string between de- 
fects opens the possibility of detecting these objects at 
low temperatures and without driving currents, simplify- 
ing the relevant experimental conditions. 




FIG. 4. Image of depicting a string connecting two de- 
fects in an experimentally realizable situation. 



Summary- In the context of the Ginzburg-Landau the- 
ory, we have given a detailed treatment of the microscopic 
aspects of pinning phenomena in nanoscale periodic ar- 
rays. Strikingly, we see an apparent first order transition 
involving the creation of a string connecting the vortices 
to the defect, and an accompanying abrupt transforma- 
tion of the supercurrent and the magnetic field it gener- 
ates. The string configuration leads to a region of lin- 
ear pinning. Absorption of the vortex by the antidot is 
marked by a jump in the pinning potential. Various as- 
pects of the pinning potential should be observable in 
experiment. 
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